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ATTENUATION IN A SHOCK TUBE DUE TO UNSTEADY-BOUNDARY-LAYER ACTION 1 


By Harold Mirels 


SUMMARY 

A method is 'presented for obtaining the attenuation of a 
shock wave in a shock tube due to the unsteady boundary layer 
along the shock-tube walls. It is assumed that the boundary 
layer is thin relative to the tube diameter and induces one- 
dimensional longitudinal pressure waves whose strength is 
proportional to the vertical velocity at the edge of the boundary 
layer. The contributions of the various regions in a shock 
tube to shock attenuation are indicated. 

The method is shown to be in reasonably good agreement with 
existing experimental data. 

INTRODUCTION 

A shock tube consists of a fluid at high pressure (region 4 
of fig. 1 (a)) separated by a diaphragm from a fluid at low 
pressure (region 1). When the diaphragm bursts, a shock 
wave propagates into region 1 while an expansion wave 
propagates into region 4. A time-distance plot of these 
waves under ideal conditions is indicated in figure 1 (b). 
Regions 2 and 3 have the same velocity and pressure but 
have different temperatures. The interface between regions 
2 and 3 is referred to as a contact surface. The analysis of 
the flow for perfect fluids is straight-forward (see, for ex- 
ample, ref. 1). In an actual shock tube, however, viscosity 
and heat conduction cannot be ignored. These lead to a 
boundary layer along the walls of the shock tube, as indi- 
cated in figure 1 (c). The boundary layer introduces non- 
uniformities into the shock tube. Analytical studies of this 
boundary layer are presented in references 2 to 6. One of 
the important consequences of the wall boundary layer is 
that it generates weak pressure waves which catch up with 
and attenuate the shock wave propagating into region 1. 
This attenuation has been studied experimentally and 
analytically in the work of references 1, 4, 5, and 6, and is 
the subject of the present report. It is assumed that the 
boundary layer is thin relative to the shock- tube diameter. 
This is a practical restriction, since most shock tubes are 
designed so that the core of potential flow is relatively 
uniform in order to permit aerodynamic tests. 

A few remarks concerning previous shock-wave-attenuation 
analyses are appropriate. In reference 4, the coordinate 
system is defined so that the shock wave is stationary. The 
flow between the shock wave and the contact surface is 
considered as a one-dimensional steady flow so that at each 


instant the mass flow through the shock wave equals the 
mass flow at the contact surface. If the mass flow at the 
contact surface is known at each instant, the correspond- 
ing shock strength can be found. The mass flow at the 
contact surface is determined from the local boundary-layer 
displacement thickness and free-stream conditions corre- 
sponding to an unattenuated shock. However, it can be 
shown that free-stream conditions do not remain constant 
at the contact surface (because of perturbations induced by 
the boundary layer). Moreover, the method of reference 4 



^Boundary layer 

(c) 


(a) Shock tube before diaphragm burst. 

(b) Wave diagram for perfect fluid, 
(c) Flow in shock tube with real fluid. 

Figure 1. — Shock-tube phenomena. 


1 Supersedes NACA TN 3278, “Attenuation in a Shock Tube Due to Unsteady-Boundary-Layer Action,” by Harold Mirels, 1956. 
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does not take into account the existence of weak pressure 
waves between the contact surface and the shock wave. 
Since it is precisely these pressure waves which are responsi- 
ble for shock attenuation, the method of reference 4 cannot 
be expected to give accurate quantitative results. It does, 
however, indicate some of the important parameters in- 
volved. Reference 1 uses the method of reference 4 in its 
studies of shock attenuation. At best, qualitative agree- 
ment with experiment is indicated. 

In reference 6, the flow perturbations due to the boundary 
layer are considered in detail. At each section of the shock 
tube, the velocity and temperature variations associated 
with the boundary layer are averaged across the tube to 
provide an equivalent one-dimensional flow. The wall 
shear and heat addition due to dissipation and heat transfer 
at the walls (all found from boundary-layer theory) are 
assumed to act on this equivalent one-dimensional flow, 
their action resulting in the generation of one-dimensional 
pressure waves propagating in both the upstream and down- 
stream directions. By integrating along characteristic lines, 
the attenuation of the shock propagating into region 1 is then 
found. The theoretical trends appeared to be in good 
agreement with the particular experimental results reported 
therein. 

Reference 5 considers the flow in a coordinate system 
which is stationary with respect to the shock. The unsteady 
nature of the flow between the shock' wave and the contact 
discontinuity is associated with the receding of the contact 
surface with respect to the shock wave. Weak pressure 
waves are assumed to be generated just in front of the 
receding contact surface. The magnitude of these waves 
is obtained by a one-dimensional averaging procedure 
similar to that of reference 6. These waves overtake the 
shock and result in attenuation. The method of reference 5 
ignores the contribution of region 3 and cannot be expected 
to yield good quantitative agreement with experiment. 

Of the previous reports on shock attenuation, reference 6 
appears to give the best agreement with experiment. How- 
ever, it can be shown that the method of reference 6 does not 
give a completely valid representation of the wave phe- 
nomena induced by the boundary layer along the shock-tube 
wall. The deficiency of reference 6 is mainly associated with 
its use of wall shear in the determination of the perturbation 
pressure waves generated by the wall boundary layer. It is 
well recognized that the effect of a boundary layer on its 
external flow is directly related to the vertical velocity at the 
edge of the boundary layer. For example, reference 7, which 
is concerned with the Rayleigh (impulsive-plate) problem for 
a compressible fluid, shows that the boundary layer generates 
pressure waves in the external flow which are equivalent to 
those which are produced if the wall moves normal to itself 
with a velocity equal to the vertical velocity at the edge of 
the boundary layer. The nonzero pressure gradient over a 
flat plate moving at high speeds (because of the finite dis- 
placement thickness of the boundary layer) is an equivalent 
steady-flow phenomenon (e. g., ref. 8). Thus, a proper way 


to find the waves generated by the unsteady wall boundary 
layer in a shock tube is to base the calculation on the vertical 
velocity at the edge of the boundary layer. Such an analysis 
is presented herein. The quantitative results for shock at- 
tenuation thus obtained would be expected to differ from 
those of reference 6. The comparison between the method 
of reference 6 and that of the present report is discussed more 
fully in the main body of the report and in the appendix 
titled GENERATION OF PRESSURE WAVES BY WALL 
SHEAR AND HEAT ADDITION. 

ANALYSIS 

One-dimensional flow with mass sources is treated, and the 
waves generated by these sources are derived. The waves 
generated by unsteady-boundary-layer action in a tube are 
found, assuming that the process can be considered as a one- 
dimensional unsteady flow. The application to the shock- 
tube problem is then indicated. 

GENERATION OF WAVES BY MASS SOURCES 

Consider one-dimensional uniform flow in a tube of con- 
stant cross-sectional area. Assume that weak mass sources, 
uniformly distributed across each cross section, are present 
and perturb the flow. Denote the net perturbation of a 
quantity (from the uniform -flow conditions) by A. (See 
appendix A for definition of symbols.) The equations of 
motion are 


/bAu , dA^\ bAp 

» (-w +w ^) = ^ ( Momentum ) 

bAp . bAu . bAp /r . , . 

-w+ p ~te +u ^ =m ( Contuiuit y) 

Ap—a 2 Ap (Isentropy) 


a) 


where m=m(x,t ) is the rate of mass addition per unit cross- 
sectional area per unit x. The presence of the sources 
generates waves. Let the superscripts + and — indicate 
perturbations associated with waves moving in the -\-x~ and 
— ^-directions, respectively. If |,r designates integration 
variables for x,t, the solution for the net perturbation at any 
point x,t may be expressed as 

Ap = A p + + A p~ (2a) 

A u= Au + +Au~ 


where 


pa 


(Ap + —Ap ) 



Ap+ ~2(1 +M)f-„ m (*’* a+u) d * 


(2b) 

(2c) 
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The integrations are conducted along the characteristic lines 
T=t— (x— £)/(a,-hu) and T=t—(£—x)/(a—u ) in the £,r plane. 
The upper limit on the integral for A p~ is <» or — «. de- 

pending on whether M< 1 or M> 1, respectively. Equation 
(2b) incorporates the acoustic relations Ap + =paAu + and 
Ap-=—paAu~. 


GENERATION OF WAVES BY UNSTEADY-BOUNDARY-LAYER ACTION 


Consider a tube of uniform cross section to have flexible 
walls such that a small normal velocity v can be generated 
at the walls. This is equivalent to mass entering the tube 

at the rate (|) pv d l, per unit x, where the integral is taken 

around the perimeter of the tube cross-sectional area. If the 
flow in the tube is considered as one-dimensional, 2 the equiv- 
alent source strength is 



4 pv 
: d 


( 3 ) 


where dz=4A/l is the hydraulic diameter, and v=v(x,t). The 
expression for A p can then be written as 

< 4 > 

Similar expressions can be written for A u and A p. In the 
case of waves induced by boundary-layer action, the v in 
equation (4) refers to the normal velocity at the edge of 
the boundary layer. Note that a positive v results in 
compression waves, while a negative v results in expansion 
waves. 

ATTENUATION IN A SHOCK TUBE 

Equation (4) can be applied to find the attenuation in a 
shock tube. The details of the analysis are described in 
appendixes B to F. The resulting formulas for shock 
attenuation are presented in the following section. The 
limiting solution for weak shocks and a numerical solution 
for an air-air shock tube are also indicated. 

SHOCK ATTENUATION FORMULAS 

The flow in a shock tube is assumed to consist of the ideal 
basic flow plus small perturbations due to the boundary 
layer. Ideal shock-tube flow relations are summarized in 
appendix G. The expansion wave of the ideal flow is 
assumed to have negligible thickness and to propagate into 

2 If the flow through the tube cannot be considered as one-dimensional, it is necessary to 
consider each element of tube surface as an elemental wave source of strength proportional to 
the local value of v. The net wave strength at any point in the tube is found from an integra- 
tion over the entire tube surface. 

In the problem of attenuation in a shock tube, the shock wave is considered to be uniform, 
laterally, at each value of x, but to decrease in strength with increases in x. Hence, a one- 
dimensional analysis is permissible. Actually, small lateral variations of shock strength 
exist, particularly near the walls, since the boundary layer immediately behind the shock 
exerts a three-dimensional effect on shock strength. These lateral variations are ignored 
herein. 

The shape of the shock wave was previously studied (ref. 9). It was assumed that region 2 
was infinite in extent (i. e., region 3 was neglected), and the problem was considered as a 
steady-flow problem in a coordinate system moving with the speed of the undisturbed shock. 
It was found that the shock wave assumes a parabolic shape, but that appreciable shock cur- 
vature is restricted to a region near the wall less than a boundary-layer thickness in extent. 


l egion 4 with velocity u— —a h as indicated in figure 2. Let 
point d of figure 2 represent an arbitrary point on the shock- 
wave characteristic. The problem is to find the net pressure 
perturbation behind the shock (i. e., A p 2 ,d)- This requires 
an integration of equation (4) along all the characteristic 
lines which contribute to A p 2td . The major contribution to 
A p 2 , d comes from characteristic lines bd, be, and ab. Hence, 
these are the characteristic lines considered in the present 
analysis (as is indicated in appendix B). Integrating along 
these characteristic lines permits A p 2 , d to be expressed as 


ad A p 2 , d 

F2y 2 p 2 



c n 

i 


©a d£+ 


T32 D 



1-CE 


A pz, d ) 


( 5 ) 


Equation (5) is derived in appendix F. The notation is 
defined in appendixes A and F. The term CE Ap 2iC jAp 2td 
represents the contribution to A p 2id of all the characteristic 
lines other than bd, be, and ab (appendix F). 

If the boundary layer in region 2 is wholly laminar or 
wholly turbulent, the boundary-layer solution for v 2 is known 
(appendixes D and E). Similarly, if the boundary layer in 
region 3 is wholly laminar or wholly turbulent, the solution 
for is known (appendixes D and E). For these special 
cases, the integrals of equation (5) can be readily evaluated. 
Thus, if the boundary layers in regions 2 and 3 are both 
laminar, equation (5) becomes (from eqs. (F6) and (F7) with 
n 2 —n z —l/2) 


ja 2 d A p 2td 

^ti 2 Tj 2 F 

\ ^ P2 n 

, l-{u s ja 2 )~\ A 

L ] 

+ M 2 J 


1 + 0 — — 
x 


&32 


r 1 1 — (u s /a 2 ) ~\ 

L m 2 J 

l+Af 3 +a 43 


>32)^32 \\}- CE (~) A Y 


( 6 ) 


Line 

Slope, 
d r/df 

oed 

\/u s 

ob 

\/U£ 

oa 

-1/(74 

ab 

l/(<7 3+ </ 3 ) 

be 

\-\/(02-e 2 ) 

bd 

l/((?2 +^2 ) 



Figure 2. — Characteristic lines appropriate for study of attenuation 

in shock tube. 
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where L 2 and X 3 are found from appendix D. If the bound- 
ary layers in regions 2 and 3 are both turbulent, equation 
(5) becomes (from eqs. (F6) and (F7) with n 2 — 7^—115) 


\uj \xj \V 2 J p 2 
f . 1 —(usla 2 ) 

Z>y 32 ML M 2 
a 32 lH-M 3 +a4 3 


-2.507 2 Z 2 F| 

\a 2 ; 

3/5 * 

r . i !-(•», wi 

W* 

L 1+ J 



] 


(y 22 ) 1/& L 22 


(— (t)T 


(7) 


where L 2 and Z 3 are found from appendix E. The other 
perturbation quantities directly behind the shock can be 
found from A p 2 ,alV 2 and normal shock relations. For 
example, utilizing equations (C4) and (C5) gives 


Au 2 , d _ 2y 1 M 2 s — (7i— 1) 1 Ap 2td 

u 2 4nfiM 2 M 2 — 1 p 2 

A p 2 , d ___ 2y 1 M 2 —(y 1 —l) A p 2 ,a 
p 2 7iMf[(7i-l)Mf+2] p 2 


^T 2>d __ Ap 2td A p 2t g 
T 2 p 2 p 2 


(Eq. of state) 


Note that the perturbations directly behind the shock are 
not isentropically related except for the limit as M s ap- 
proaches 1. This is due to the fact that a perturbation of 
shock strength creates an entropy perturbation. (The 
particles between the shock wave and the contact surface 
have entropy perturbations whose magnitude depends on 
the strength of the shock at the instant each particle passed 
through it. To the order of the present analysis, these 
entropy perturbations do not interact with the longitudinal 
wave system induced by the wall boundary layer.) 

Expressions similar to equations (6) and (7) can be found 
for the mixed cases n 2 = 1/2, n 2 = 1/5 and n 2 =lf5, 7& 3 =l/2. 
For other boundary-layer characteristics (such as transition 
from laminar to turbulent flow in the middle of regions 2 or 
3), it is necessary to integrate equation (5) with the corre- 
sponding v distribution. Reliable criteria for determining 
the transition points in shock-tube boundary layers have 
not yet been established. A crude tentative method for 
estimating the transition points in regions 2 and 3 is pre- 
sented in appendix H. It is pointed out therein that the 
Reynolds numbers at point b, computed separately for re- 
gions 2 and 3, can be used as an index to determine whether 
the boundary layers in regions 2 and 3 are primarily laminar 
or turbulent, respectively. The values of Re 2lb and Re Sfb , as 
defined in appendix H, can be obtained from figure 3 for the 
case where the gas in region 1 is air at 520° R. 

The present theory requires that the boundary layer be 
thin relative to the shock-tube diameter. The boundary- 
layer thickness in region 2 is defined by equations (D5) and 
(E3) for laminar and turbulent boundary layers, respectively. 
The laminar-boundary-layer thickness S 2 is taken to corre- 
spond to u/u 2 — 0.99, while the turbulent-boundary-layer 



Figure 3. — Reynolds number at point b for air in region 1 at 520° R. 
Rez.b=vnRe 2 . b (appendix H). 


thickness 8 2 is taken to be the value obtained from a Karman- 
Pohlhausen-type integral solution. Considering conditions 
along characteristic lines, the maximum boundary-layer 
thickness occurs at point b of region 2. Values of d 2tb and 
8 2<b are plotted as a function of M s in figure 4 for the case 
where the gas in region 1 is air at 520° R. 

It is previously noted, in connection with equation (4), 
that a positive v results in the generation of compression 
waves, while a negative v results in the generation of expan- 
sion waves. From the boundary-layer theory of appendixes 
D and E it can be seen that, apart from dissipation and 
heat-transfer effects, v 2 is negative and v z is positive. Thus, 
the boundary layer in region 2 induces expansion waves 
(which attenuate the shock), while the boundary layer in 
region 3 induces compression waves (which accelerate the 
shock). Dissipation and heat transfer modify these results. 
Dissipation tends to increase v in both regions 2 and 3. 
Heat transfer from the wall to the boundary layer increases v. 
Heat transfer from the boundary layer to the wall decreases v. 
In region 2 the heat transfer is from the boundary layer to 
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Nominal shock Mach number, M s 


(a) Laminar case (eq. (D5)). 



Nominal shock Mach number, M $ 

(b) Turbulent case (eq. (E3)). 

Figure 4. — Boundary-layer thickness at point 6 in region 2. Gas in 
region 1 is air at 520° R; wall is assumed a perfect conductor 

(r lt .= T 1 ). 


the wall (ref. 2) and leads to larger negative values of v 2 and 
therefore more shock attenuation. In region 3 the heat 
transfer is from, the wall to the fluid for a weak expansion 
wave and from the fluid to the wall for a strong expansion 
wave (ref. 2). Thus, for a weak expansion wave, heat 
transfer in region 3 tends to generate compression waves 
(thus accelerating the shock). For a strong expansion wave, 
the heat transfer in region 3 tends to induce expansion waves 
(which attenuate the shock). The net effect of all these 
factors is to attenuate the shock wave. The relative magni- 
tudes of the various terms in equations (6) and (7) are noted 
later in a numerical example. 

LIMITING SOLUTION FOR WEAK SHOCKS 

For M s approaching 1, equations (6) and (7) take on the 
following forms: 

Laminar case: 


(Us d a 2 d\ V2 Aj? 2 ,d __ — 8 Ti A j 7i— 1 \ 

\u 2 X V 2 ) P2 -Vx(Yi+l) V 


(for T 2 , U =T X ) 

(8a) 


= 7 === (for T,. „= r 2 , r ) (8b) 

Vtt(7i + 1) 


Turbulent case: 



a 2 d\ A ‘Ap 2 ,d_ 


— 0.11 57? 
3 [2(Ti+1)]^ 


(9a) 


-0.1157! 
[2(7i + 1)]^ 


(for T 2 ' W =T 2 ,t) (9b) 


The condition T 2tW ~Tx (eqs. (8a) and (9a)) corresponds to 
the case where the shock-tube wall is a perfect conductor 


(assuming that the fluid in region 1 is in thermal equilibrium 
with the wall). Note that taking T 2tW equal to T x means 
that the wall remains at its original temperature and that 
heat is transferred from the fluid (of region 2) to the wall. 
The condition T 2iW =T 2 , t (eqs. (8b) and (9b)) corresponds to 
the case wherein the wall is a perfect insulator or has a very 
low heat capacity. When the shock-tube wall is a metal, the 
assumption T 2>w =Ti should give ver 3 r accurate results, as 
discussed in reference 2 (particularly since weak shocks are 
now being considered). Results for both T 2tW =Tx and 
T 2tW =T 2>r are given so as to define the upper and lower 
bounds of the heat-transfer effect on shock attenuation. 
The attenuation is greater, by a factor of approximately 
7i , when there is heat transfer as compared with the insulated- 
wall case. (The factor is exactly 7i for the turbulent 
case and approximately 7i for the laminar case, the latter 
depending on the value of <ii). 

AIR-AIR SHOCK TUBE WITH T 4 =Ti=520° R 

Equations (6) and (7) were evaluated numerically for an 
air-air shock tube with T 4 =Ti=520° R. It was assumed 
that cr=0.70 and 7—1.4. The results for the laminar and 
turbulent cases are summarized in figure 5. Both the 



Nominal shock Mach number, M s 


(a) Laminar case (eq. (6)). 



Nominal shock Mach number, M s 

(b) Turbulent case (eq. (7)). 

Figure 5. — Attenuation in air-air shock tube. T l ~T A = 520° R; 
Prandtl number <r, 0.70; ratio of specific heats y, 1.4. 
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Figure 6. — Percent contribution to shock attenuation of characteristic lines of figure 2. Air-air shock tube; T\—Ti— T 2 ,w~ T 3 >U) = 520° R; 

Prandtl number er, 0.70; ratio of specific heats y , 1.4. 


insulated-wall case (T 2w — T 2 y,T$ tW = T s>r ) and the case where 
the wall is a perfect conductor {T 2tW =T ZtW =T 1 =T^) are 
noted therein. When the shock-tube wall is a metal, the 
assumption that the wall is a perfect conductor should give 
accurate results except possibly for very strong waves (e. g., 
ref. 2). 

The relative contribution to the net attenuation of the 
various terms in equations (6) and (7) is indicated in figure 6 
for the T 2 , w =Tz,w=T 1 =T 4; case. In particular, the percent 
contribution of the integrations along the characteristic lines 
bd , ab, and be and of the reflected wave at point c are indi- 
cated therein. (The reflected wave at point c represents the 
contribution of all the characteristic lines in figure 2 other 
than lines bd, ab, and be, as is mentioned in appendix F.) 
For weak shocks, the major contribution to shock attenua- 


tion comes from the integration along line bd. With increas- 
ing M s , the contribution of the characteristic line be increases 
gradually to a value of about 30 percent at M s = 6.0. The 
integration along line ab leads to compression waves (which 
tend to accelerate the shock) and therefore is negative in 
figure 6. Its value decreases to about —20 percent at 
M s ^2.0 for the laminar case and to about —35 percent at 
M s &i 2.75 for the turbulent case, and then increases with 
increases in M s . (The influence of characteristic line ab 
may be somewhat overestimated in the present analysis 
because of the assumption of an expansion wave of zero 
thickness.) The contribution of the reflected wave at point c 
varies from a value of zero at M s — 1 to about —5 percent at 
M s — 6. The neglect of the latter contribution appears 
reasonable for the range of M s considered herein. 
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RESULTS AND DISCUSSION 

Equations (6) and (7) define the attenuation in a shock 
tube when the boundary layer is wholly laminar or turbulent, 
respectively. These equations are now compared with the 
experimental and theoretical results of reference 6. 

COMPARISON WITH EXPERIMENTS OF REFERENCE C 

Measurements of shock attenuation were obtained in the 
investigation of reference 6 by using a high-pressure shock 
tube having a %- by %-foot rectangular cross section. Air at 
room temperature was used in regions 1 and 4. The air 
in region 1 was maintained at atmospheric pressure. Four 


sets of runs, corresponding to p 41 =4.061, 5.764, 7.455, and 
17.915, are reported. The results of these runs are summa- 
rized in figures 7 (a) to (d), respectively. Some of the theo- 
retical curves of reference 6 are included in these figures. 
The theoretical predictions of equations (6) and (7) are also 
indicated in figure 7. The latter were found from figure 3 
with Pi= 2117, and by assuming the wall to be a perfect 
conductor (since the shock tube had metal walls except for a 
pair of schlieren glass inserts) . 

The Reynolds number at the contact surface for these 
tests is given directly by figure 3. The boundary-layer 
thickness at point b of region 2 for the test conditions of 



(a) Pressure ratio p t /p h 4.061; shock Mach number relative to wall (b) Pressure ratio p*/p u 5.764; shock Mach number relative to wall 

M„ 1-344. m s , 1.442. 
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(c) Pressure ratio pt/pi, 7.455; shock Mach number relative to wall (d) Pressure ratio p*/pi, 17.915; shock Mach number relative to wall 


M s , 1.518. 


M s , 1.792. 


Figure 7. — Pressure behind shock wave as function of distance from diaphragm. Air-air shock tube. Ti=Ti= 520° R; hydraulic 

diameter d, 1 /7. 
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Figure 8. — Relative boundary-layer thickness at point b of region 2 
for experimental conditions of reference 6 (air in region 1 ; turbulent 
boundary layer; T l =T 2tU) =520° R; h, 1/8 ft). 


reference 6 is plotted in figure 8. It is assumed in figure 8 
that the boundary layer in region 2 is wholly turbulent and 
that the wall is a perfect conductor. The boundary-layer 
thickness is presented in the form of a relative thickness 
2 82 , 1 ,! h where h is the smallest shock- tube dimension normal 
to the flow and equals % for the shock tube of reference 6. 
When 2d 2 ,b/h is small compared with 1, the assumption of a 
thin boundary layer is valid, and the theory of the present 
report is applicable. It may be seen from figure 8 that 
25 2 ,&/^ 0.20 for r=S12 for the test conditions of reference 6. 
Since 25 2 , blh tends to overestimate the effective boundary- 
layer thickness (because of the asymptotic manner in which 
the velocities approach free-stream conditions at the edge 
of the boundary layer), 3 * figure 8 indicates that the theory 
of the present report is applicable for comparison with the 
experiments of reference 6. 

In figure 7 (a), the first data point (as; £*3') is close to the 
theoretical value for a wholly laminar boundary layer. The 
other points fall somewhat above the theoretical curve for a 
wholly turbulent boundary layer. The Reynolds number 
per foot at the contact surface for this case is Re 2yb fx=0Ay. 
10 6 (fig. 3). Using I&?=0.5X10 6 as a rough indication of 
the transition Reynolds number, as discussed in appendix H, 
it is reasonable to have the first data point near the laminar 
curve and the other points near the turbulent curve. 

The experimental results for M x = 1.442 are given in figure 
7 (b). The Reynolds number at the contact surface is 


3 A properly computed displacement thickness would probably give a better estimate of 

the effective boundary-layer thickness. 


Re 2 tb (x= 0.8 X10 6 . The data agree quite well with the 
values for a wholly turbulent boundary layer. There is a 
slight tendency for the points to lie above the theoretical 
curve, which might be attributed to the short length of 
laminar boundary layer directly behind the shock wave. 

The experimental results for M s = 1 .518 ( Re 2 yb jx — 1 . 1 X 10 6 ) 
are given in figure 7 (c) and are in excellent agreement 
with the turbulent-boundary-layer theory of the present 
report. 

In figure 7 (d), the experimental results {M s — 1.792, 
Re 2 >6 /r=2.8X10 6 ) fall somewhat below the theoretical pre- 
dictions of the turbulent-boundary-layer theory of the 
present report. 

In general, the theory seems to agree reasonably well with 
experiment for the range of data considered in figure 7. 
Figure 7 (d) indicates the poorest correlation and suggests 
that the present theory may underestimate the attenuation 
corresponding to large values of M s . 

COMPARISON WITH THEORY OF REFERENCE 6 

The present analysis assumes a relatively thin boundary 
layer, and the calculations are based on the vertical velocity 
at the edge of the boundary layer. The numerical value of v 
depends on terms which are related to the boundary-layer 
velocity profile, dissipation, and heat transfer. See, for 
example, equations (D3a) and (D3c). These terms cor- 
respond to the use of wall shear, dissipation, and heat trans- 
fer in reference 6. Dissipation and heat transfer play similar 
roles (qualitatively) in the present analysis and in reference 
6. However, it can be shown that the velocity-profile term 
in the equation for v has a sign opposite to that of the wall- 
shear term in reference 6 for characteristic lines ab and be 
(e. g., appendix I). Therefore, the integrations along lines 
ab and be, in reference 6, overestimate and underestimate, 
respectively, the contributions of these characteristic lines 
to shock attenuation. If the same boundary-layer theory 
is used, the attenuation calculations of the present report 
and of reference 6 should be in quantitative agreement 
only in the limiting case M s ^ 1 (for which the contributions 
of characteristic lines ab and be are negligible) and for those 
values of M s where the errors due to characteristic lines 
ab and be tend to compensate. 

Since different boundary-layer theories were used in 
reference 6 and herein, it is not directly possible to separate 
the discrepancies between the two methods due to the re- 
spective boundary-layer theories from the discrepancies due 
to basing the attenuation calculations on wall-shear, dissipa- 
tion, and heat-transfer terms rather than on v. An estimate 
of the latter discrepancy can be obtained by reversing the 
sign of the first term in the equations for L 2 (eq. (D3c)) 
and Zr 3 (eq. (D4c)) for the integrations of v along character- 
istic lines be and ab. Such a procedure shows that the theory 
of the present report and that of reference 6 agree at M s —\. 
With increasing M s , reference 6 first overestimates the at- 
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tenuation because of the increasing importance of character- 
istic line ab. At M s =0(1.5), reference 6 overestimates at- 
tenuation by about 10 to 15 percent. With further increases 
in M s , the errors in characteristic line be become important 
and tend to compensate for the errors in line ab so that the 
discrepancy between the two methods decreases. At 
M s = 0(2), the two methods are again in approximate agree- 
ment. As M s increases further, line ab becomes relatively 
less important compared with be, and the theory of reference 
6 underestimates the attenuation. At M s = 6, the method 
of reference 6 appears to underestimate the attenuation by 
about 50 percent. The above figures are only approximate 
because of the manner in which they were obtained, but they 
indicate the proper trend with M s . 

From figure 5 (b) it can be seen that the attenuation theory 
of reference 6 (for turbulent boundary layers) is in good 
agreement with that of the present report for 1<M S <2. 
The agreement is somewhat better than that to be expected 
from the discussion of the previous paragraph. Hence, 
either theory could be used to correlate the experimental data 
of figure 7. For M $ ^> 2, reference 6 considerably underesti- 
mates the attenuation. The laminar-boundary-layer theory 
of reference 6 differs by a factor of about 2 from that of the 
present report; therefore, there is a large discrepancy between 
the theoretical laminar-boundary-layer curves in figure 7 (e) . 

CONCLUDING REMARKS 

A method is presented for computing the attenuation of a 
shock wave due to unsteady-boundary-layer action. The 
various assumptions involved in the analysis are summarized 
in this section, since these define the limitations of the 
method and suggest possible fields for improvements: 

1. Small perturbations: The equations of motion were 
linearized, assuming the potential flow external to the wall 
boundary layer undergoes only small perturbations. For 
long shock tubes with large amounts of attenuation, it might 
be advisable to employ a characteristic method. 

2. Thin boundary layer: The assumption of a thin bound- 
ary layer (relative to shock-tube diameter) is consistent with 
the assumption of small perturbations of the potential flow. 
If the viscous effects span the entire tube cross section (i. e., 
long tubes), it may be advisable to base the shock-attenuation 
theory on wall shear. 

3. One-dimensional longitudinal waves: It was previously 
pointed out that the perturbations in a shock tube should be 
computed by assuming each element of wall surface area to 
be an elemental acoustic source of strength proportional to v. 
This gives rise to a complex wave pattern involving both 
longitudinal and transverse waves. For the purposes of the 
present analysis it was assumed that the longitudinal waves 
are of primary interest and that these can be computed on 
the basis of a simplified one-dimensional theory. This 
assumption is accurate when an observer is relatively far 


from the sources (since the details of the source distribution 
around a perimeter then become relatively unimportant) but 
introduces errors when the observer is near the sources. It 
would seem that the assumption of one-dimensional longi- 
tudinal waves is accurate for weak shocks but introduces 
errors for stronger shocks, particularly when the flow relative 
to the wall is supersonic. 

4. Boundary-layer theory: The boundary-layer theory of 
reference 2 was used. It may be assumed that the laminar- 
boundary-layer solution is reliable except possibly for very 
strong shock waves. In the latter case, very large tempera- 
ture gradients exist normal to the wall, and it may be advis- 
able to choose a different reference temperature from that 
used herein. Also, for strong shock waves, dissociation might 
occur, which would require changes in the laminar-boundary- 
layer theory. The turbulent-boundary-layer theory of 
reference 2 requires experimental verification for even the 
weak-shock case. However, the good agreement between 
the attenuation calculation based on the turbulent bound- 
ary layer of reference 2 and the experiments of reference 6 
suggests that the turbulent-boundary-layer theory of refer- 
ence 2 gives reasonable results, at least for the weak-shock 
case. 

5. Expansion wave of zero thickness: In order to simplify 
the problem of determining the boundary layer behind an 
expansion wave, it was assumed in reference 2 that the 
expansion wave is of negligible thickness (i. e., “expansion 
shock”)- This assumption is valid for weak expansion 
waves but is in error for strong expansion waves. The 
boundary-layer solution for region 3 might be improved for 
the strong-wave case by considering the finite thickness of 
the expansion wave. However, the contribution to shock 
attenuation of region 3 becomes small for the strong-wave 
case, so that an improved boundary-layer solution for region 
3 may not significantly affect the attenuation calculation. 

In addition to the preceding discussion, the following 
extensions of the present report might be pursued: 

1. The details for obtaining the perturbations at a fixed 
point in the shock tube, as opposed to finding the shock 
attenuation as is done herein, might be treated. This would 
be useful for further correlating the theory with experiment 
and for determining conditions at an aerodynamic model 
(when the shock tube is used as a wind tunnel). 

2. Formulas equivalent to equations (6) and (7) might 
be derived for the case where the boundary layer is partly 
laminar and partly turbulent. 

3. The boundary-layer theory as outlined in appendixes 
D and E might be simplified, particularly for the turbulent 
case, so as to give reasonably accurate results without too 
tedious a boundary-layer calculation. 

Lewis Flight Propulsion Laboratory 

National Advisory Committee for Aeronautics 
Cleveland, Ohio, April 24, 1956 
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APPENDIX A 

SYMBOLS 


The following symbols are used in this report: 

A cross-sectional area of tube 

a speed of sound 

C eq. (C3) 

c P specific beat at constant pressure 

c v specific heat at constant volume 

D eq. (C3) 

d hydraulic diameter, AAfl 

E eq. (C7) 

F eq. (C7) 

h smallest shock-tube dimension normal to flow 

k thermal conductivity 

L 2 ,L 3 appendixes D and E 

l perimeter of shock-tube cross section 

M Mach number of flow relative to wall 

Mg shock Mach number relative to wall 

M 2 ,M 3 Mach number of flow in regions 2 and 3, relative 
to wall 

m rate of mass addition per unit cross-sectional area 

per unit x 

n 2l n 3 appendixes D and E 

p pressure, lb/sq ft 

Re Reynolds number (appendix H) 

Re 2 ,b Reynolds number at point 6 of region 2, 

v WSJ( 1+M2 -S) M *} appendix H 

T temperature, °R 

T r temperature of insulated wall 

t time 

u velocity of flow relative to wall 


u 9 velocity of shock wave relative to wall 

u 2 ,u 3 velocity in regions 2 and 3 relative to wall 

v vertical velocity (positive when directed into tube) 

at edge of boundary layer 
* longitudinal distance 

7 ratio of specific heats 

A perturbation quantity (Ap= perturbation of p, etc.) 

8o laminar-boundary-layer thickness in region 2 (eq. 

(D5) ) 

<5 2 turbulent-boundary-layer thickness in region 2 (eq. 

(E3)) 

n coefficient of viscosity 

v kinematic viscosity 

£ integration variable representing x 

p mass density 

<r Prandtl number 

r integration variable representing t 

Subscripts : 

1,2, 3, 4 regions of shock tube (fig. 1) 

a,b,c,d,e points on characteristic lines (fig. 2) 
m evaluated at mean temperature of boundary 

layer (appendixes D and E) 
w evaluated at wall 

Superscripts : 

+ associated with wave moving do wnstream 
(+£-direction) 

— associated with wave moving upstream 
(— ^-direction) 

Special notation : 

a 32 s a 3 /a 2 ; Pi 2 =Pi/P 2 , etc. (Two successive integer sub- 
scripts, not separated by a comma, represent a ratio.) 


APPENDIX B 


CHARACTERISTIC-LINE GEOMETRY AND SOME INTEGRALS 


The equations of the characteristic lines considered in the 
present report are summarized. Some important integrals 
are indicated. It is assumed that the characteristic lines are 
straight (corresponding to the ideal-flow case) and that the 
expansion fan can be considered as an “expansion shock” 
(i. e., expansion wave of zero thickness) moving with velocity 
a 4 into region 4 (following sketch). The point d with coor- 
dinates x,t is assumed to be on the shock-wave characteristic. 


t,r 



The equations of the characteristic lines are 
Line bd : 

£=x— {(h+U 2 )(t— r) 

x-a 


Line ab : 


Line be: 


r=t- 


U S T~%=\ 


Q'2~\"'U'2 

|~I -\-M 2 — ( u s fa 2 ) 


1-\~M 2 

£— &“-(as+tt 8 ) (r 6 — r) 


-1 (*-«) 


T = T b ~ 




■a 8 +u 3 

, ^ /lT-M 3 4-a.43\ y \ 

aiT+l= ( i'qnf 3 


£=£&+ {a*— u 2 ) (r b ~ r) 


T — Tl) - 


, ri -m 

U s T—£=^ 


(Jb 2 U 2 
1 — M 2 A (Usldz) 

M 2 


■]( fc -0 


(Bl) 
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In terms of x,t, the coordinates of points a, b, and c are 

Point a ; fe- 1 +M 3 +a a J u,\ 


Point b : 


Point c : 


. r 

a L 1+^3 + ® 43 J 

f,=a[l+Af,- («,/<*,)] ^ (B2) 

t 6= t [ 1 + M 2 — (u s ja2) ] 

r 1 -\-M-2 — {Usjfy ) ”] 

ic ~ x Ll-M*+(u s /at)J 

r 1 -{-M2 — ( u s 1 <12) ~] 

Tc \_l-~M 2 ~\- (u s la 2 ) J 

Some important line integrals are 




1 -j— 7VdT 3 -j— CX 43 \ 1 — ns) 


(S3) 


Weak pressure waves are assumed to overtake a contact 
surface or a shock wave. The strength of the reflected and 
transmitted waves is found. 

CONTACT SURFACE 

Consider the contact surface which separates region 3 from 
region 2 (following sketch). At a certain instant, known 


APPENDIX C 

INTERACTION OF PRESSURE WAVES WITH INTERFACES 

where 



incident waves Ap^ b and Ap^b intersect the surface. The 
problem is to find the final waves Ap 3 ~ 6 and Ap} b . From 
isentropic flow relations, 

Ap + =pa Au + 

Y (CD 

A p =—pa Au J 

The boundary conditions across the contact surface are 

Aul* + A ui b = A ut.b 4- Ant, b 

Ap£b+Api'i,=Ap£ b +Apl b 

From equations (Cl) and (C2), the expressions for Apl & and 
A pi b are 

Ap£ b = CApl b + Z? A pi b 
A » =7320-23# Apib~CAptb 


G~ (732®23~ 1)/(Y32®23H“ 1) 

D~ 2/(732^23 + 1) 

The quantities C and D are reflection and transmission co- 
efficients, respectively. 

SHOCK WAVE 

Consider a shock wave moving with velocity u 3 into a 
stationary fluid (region 1). The region behind the shock is 
designated as region 2. Conditions in regions 1 and 2 can be 
found from shock-wave theory. In particular, with y 1 =y 2 , 

2YiA^— (7i— r> 


2fci= 

u 2 


7j + l 


(*-£) 


dl 7 l + 1 

(7i+l )M 2 8 


(04) 


(C2) 


(03) 


p 21 ”'(7i-l)i^f+2 

At a certain instant, a known wave Ap£ d intersects the 
shock (following sketch). 
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The problem is to determine the net perturbation of the 
shock wave. From equations (C4), 




A u 2(d 


ai 7i+l 

Ap 2 ,d 


( 1+ A) 


A Ms, rf 
AM,. d 


P2 (Yt-DMf+2 Ms 


(C5) 


The net pressure and velocity perturbations in region 2 are 
then related by 


Ml 

A^2.d=2pia 1 A u 2td 


(C6) 


With Ap 2 ,d=A^ d ~YAp^ id and Att^sAttfi + Au,,* using- 
equations (Cl) gives 

A^ d =J£A;p 2 %y 

Ap2,d = F'Ap2,d J 


(C7) 


where 


*-( 2pu “” rSf-O/C 2 ^ 12 ™ +1 ) 

^-( 4 ^-rS)/( 2 ^OT +1 ) 

The quantity ^ is a reflection coefiicient. The resulting per- 
turbation of the shock Mach number is 


AM s , tf = 7i + l p Ay^ 

Ms 47! Ml p 2 


(08) 


For M s = 1 , 


A^>_ 


Ap 2 .d 


=0. For M s infinitely large, 

a 


Apt 


which becomes 


AP2.C 


Apia 


‘ 2 W +1 

0.14 for 7 i= 1.4, indicating that the 


reflected wave is relatively small and is opposite in sign to 
the incident wave. 


APPENDIX D 

LAMINAR BOUNDARY LAYER BEHIND WAVE 


The laminar boundary layer behind a shock or thin 
expansion wave is analyzed in reference 2 . Some of the 
results are summarized herein. It is assumed that the wall 
temperature behind the wave is constant and that <r and c v 
are independent of temperature. The fluid properties p 
and k are referenced to a mean temperature as discussed in 
appendix C of reference 2 . The notation of reference 2 is 
compared with the notation of the present report in the 
following equations. The left side of each equation repre- 
sents the notation of reference 2 , while the right side is in 
the notation of the present report. 

Region 2 


Region 3 


V/ w — V/ S 

u w 

1 

u e =u s —u 2 

U e 1— (u 2 /u s ) 

V e =V 2 

^£_1 = 

u 2 ju s 

X—U s T—£ 

u e 

1 — (U 2 /Ug) 

U-id — H 4 

u w 

1 

u e = CU+Uz 

u e l + (u 3 /ai) 

v e =v 3 

Uw__ i = 

—Us/di 

x—a^T-j- £ 

u e 

1 ~f" (^ 3 /^ 4 ) 


(Dl) 


(D2) 


VERTICAL VELOCITY AT EDGE OF BOUNDARY LAYER 

The expression for the vertical velocity at the edge of the 
boundary layer can now be written as follows: 


Region 2 


where 


_ _ T r ui/u s "i 1 wg ( V 2 X 2 

Vz 2 Ll — (u 2 /u s )J \u t r—z) 


- (u 2 /u s ) 
n 2 =Y2 


(D3a) 

(D3b) 


JpUlim \_Jid _ 

\ 2 P 2 T 2 , m • | 7?->oo L_ ju s - _ 


7 2 


!m '(s)( 1 'S)J."' ,j ' + 


(D3c) 


fa-,) (i-?)"| 

S J- 113< Va 


lim 

>H> 0 ° 

J. 


/; 


r 2 d^= 1.569 
s 2 d u =1.134^Q 


/ .. l-(g 

V 1.933— 


(WaJUt) 


— (u 2 /u s ) 


(u 2 /u s ) 
(0-2) 


1 1 — (u 2 ju s ) 

(D3d) 

2 . 022 — (u 2 lu s ) 

“ —OMbw/u, 


- („,) !-<»>/».) 

(D3e) 

-0.504-0.47 (tts/tt.) 



022— (u 2 fus) 

0 37—0.38 (<«/«,) 

r 2 ( 0 ) = (<r 2 ) 
r i\, , , 7 2 - 1 


l- (Will,) 


^=1- 


M|r 2 ( 0 ) 


P2 


T< 


2 ,m_ 


P 2 , m 


°- 5 (% ?+1 ) +0 - 22 (% ?_1 ) 


(D3f) 

(D3g) 

(D3h) 

(D3i) 


The above quantities can be evaluated if the wall surface 
temperature T 2yW is known. The value of T 2 , w depends on 
the heat capacity and conduction properties of the wall and 
the wall thickness. A method for evaluating T 2 , w is pre- 
sented in reference 2. From the discussion in reference 2 , 
it may be concluded that in most cases T 2tW approximately 
equals T x when the wall is a metal. 
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Region 3 


Vz — La ^ 


where 


Uj/(h \ 
1 +Uz/aJ 

nz=y 2 


}- ns ( Va Y 

\a 4 T + £/ 


** +{lim f <»-■” (1+^ 1+ 

2 Hz Tz,?n V.ij->coL Uzldi J 

-M 3 a 4 3^1+^J^ T 3 d*?+ 


Ts— 3 


1 + (dajad ( T-z, i 

Uzldi 


\Tz 


%)Jo’ S3d ’’} 


r 3 d? 7 = 1.686 


vi 


i+(W g *) 


234+ (Uz/ad 


s; 

f' , '-'WiaM w 


326+ (Uzjdi) 

0.374-0.50 (tt 3 /a 4 ) 

r 3 ( 0 ) = (<r 3 ) 1+( “ 3/ ° 4) 

1 


Tz,r 1 I T 3 

rp — IT 


M 3 r 3 ( 0 ) 


Tz,m 

-0 5 w + 

1^4-0 22 — 1^) 

Tz “ 

-0.5 \T 3 j 

V+ V Tz v 


(D4a) 


(D4b) 


(D4c) 


Iimr ( ^~ /) - ( - 1 7 +M3/a<) 1^1.217J 9 . (D4d) 

*->«+ ^ 3 /«4 J v 2-326 +(^ 3 /a 4 ) 


326+ (^ 3 /^ 4 ) 

0.22U3/a< 

(<r 3 ) 1 +^ tf4 ) (D4e) 


— 0.50— 0.36(tt3/g<) 
14-(«3/fl4) 


(D4f) 


(D4g) 


(D4h) 


(D4i) 


Equations (D4) can be evaluated if T z>w is known. If the 
wall is initially at temperature T 4 , then T z , w = T . 4 is generally 
a good estimate for T 3iW , as is discussed in reference 2 . 

BOUNDARY-LAYER THICKNESS IN REGION 2 

In the body of the report it is assumed that the boundary 
layer is thin relative to the shock-tube diameter. It is there- 
fore of interest to present an expression defining boundary- 
layer thickness. Region 2 is of particular interest since the 
boundary layer in region 2 is generally larger than that in 
region 3 . Let S 2 represent the laminar-boundary-layer thick- 

u 

ness in region 2, defined so as to correspond to —=0.99. 

u 2 

From equation (18) of reference 3, 


5 2 = 


V 2 l fJ- 2 , m T 2 jv 2 {u s r— g) p .89 Vl— u 2 /u s t 

Y ^T 2 ,m \ u s Lvi*413™ u 2 Iu s " r 

u s 

d »J (D5) 


Ta— 3 


(Eq. (D 5 ) uses a mean reference temperature and an inter- 
polation formula for rjs (defined in ref. 3) and takes the upper 
limit on the integrals to be 00 rather than 775 .) To evaluate 
8 2 along the characteristic line bd of figure 2 , take u s t—£ equal 

l++f 2 UsfaL^ (x—Q' The value of b 2 at point b is found 

1 -\-M 2 — u s ja 2 \ ( , u 2 \ _ 

1+M 2 A uJ X - 


to 


(-■ 


1+M 2 

by replacing tt s r— £ 


by Q 


APPENDIX E 

TURBULENT BOUNDARY LAYER BEHIND WAVE 


The turbulent boundary layer behind a shock or thin ex- 
pansion wave is also studied in reference 2 for the case of 
constant wall surface temperature. The correspondence 
between the notation of reference 2 and that of the present 
report was previously noted in equations (Dl) and (D 2 ). 
The turbulent-boundary-layer solution of reference 2 (re- 
lating to v ) is now summarized. It is assumed that the wall 
surface temperature is essentially constant. 

VERTICAL VELOCITY AT EDGE OF BOUNDARY LAYER 



Region 2 
where 


v 2 — — L 2 


ulju s ~ l ~ n 2 / v 2 \ n 2 
_1—{u 2 /u 8 )J \U s T—%/ 


7 l 2 — 1/5 


^2 

k 


L 2 =0.0460 


8 2 [(m 2 1 / ii s ') 1] 

k(u 2 /u s ) 


— T 2 u 2 

( i u , 2 ju s ) u 2 L\T 2 , w U s 


f b 2 (u 2 lu s ) \ 4/5 

OO-’f-. 7 ")- 


7c 2 



(L 2 , 7“ 



( u *! u s~L) _ t 2 /t 2<w r / 1 „ t 2 j \ 

82 (%Mr) l—(u 2 /u s ) L\ ‘T %9 2,7 J 


(Ela) 

(Elb) 

(Elc) 


(Eld) 



-(%-')(£) 

(Eli) 



(Elj) 

T 2 ,m 

t 2 " 

=0 - 5 (% ?+1 ) + °- 22 (% r “ 1 ) 

(Elk) 


The conditions under which T 2iW essentially equals Ti can be 
established by the methods of reference 2. The assumption 
of T 2 iW =T 1 appears reasonable for most cases. 

The integrals I 2j7 and I 2 t 7 —I 2t3 can be evaluated using 
table I. Reciprocals are tabulated therein so that linear 
interpolation is accurate except for b 2 near —1 or c 2 /(& 2 + 1) 
near 1. The integrals can be evaluated analytically from 
the following expressions: 



14 


REPORT 1333 — NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 


-L,7 = 


V&2 + 4c 2 


[(4, 7 )^ 2 " {h, 7 ) x 2 ] 


-L,7 1%, 8 = 

where 


v&!+ 4c< 


[(X 2 1) (/ 2f 7 ) X 2 (02 1) (-L, 7)/3 2 ] 

^2= 2 “ (&2+V^H"4c 2 ) 
ft=2^ 0 2 -V*i+4c 2 ) 


(Ell) 

(Elm) 

(Eln) 

(Elo) 

(Elp) 


It can be shown that /3 2 <0, X 2 >1. Equations (Elo) and 
(Elp) give accurate results for — 1.6<0 2 <(O and 1<X 2 <2, 
respectively. For values of /? 2 and X 2 outside these ranges, 
the right sides of equations (Elo) and (Elp) become the 
difference of two nearly equal numbers and the following 
expansions are useful: 

*^-s^[l+S(ISi!(l= S )'] (Elo) 

(/, 7 )' s= 8( 1 -w[ 1 ~ 8 Si(7+mK8+m)(^) ] (Elr) 


Region 3 


< E2 ‘> 


where 


n z =% 


(E2b) 


L 3= 0.0460 /« (E2c) 

5 3 (« 3 K> \ ^[l + ^ayMj 

5 3 l-\-(uzM _a 4 \~ f7 ^ T z (T T ^ 

~ ~7T I * ^3 Tri U3.7 —J 3,8v , ~" 

5 3 Uz/CLi Vj Z L -*3, to 


03 1 ~f~ (^ 3 /^ 4 ) T z fT Z ' 


u z /a,4 


7 T \ 
* ™ -<3,7 1 


l 3, to 


I+m^K 1 

V/t 3 / 


xn ~L 


z N dz 


T : 


T z 


l-\-b z z—c z z 2 
63 = (T 3 , r /r 3 , «,) 1 

C 3=[(r 3 .r/r 3 )~i](r 3 /r3 fW ) 

^=l+^M§(<r 3 )* 

=°- 5 (% ?+1 ) +0 - 22 (% _1 ) 


(E2e) 

(E2f) 

(E2g) 

(E2h) 

(E2i) 

(E2j) 

(E2k) 


The conditions under which T 3tW essentially equals T 4 can 
be established by the methods of reference 2. The assump- 
tion T ZtW =T A appears reasonable for most cases. The 
integrals / 3i7 and / 3t7 — / 3> s can be evaluated using table I or 
equations (Ell) to (Elr) (with the subscript 2 replaced by 
the subscript 3 therein). 

BOUNDARY-LAYER THICKNESS IN REGION 2 

The symbol <$ 2 represents the turbulent-boundary-layer 
thickness as obtained from an integral (Karmdn- 
Pohlhausen) type of boundary-layer solution. From refer- 
ence 2, it can be shown that 


g 2 = 0.0574 4 




<P2 


u 2 /u s 


§2,(U2_ 1 \ 
8 2 \Us J 


(1 




(E3) 


The value of 8 2 along characteristic line bd or at point b may 
be obtained by the substitutions indicated after equation 
(D5). 


APPENDIX F 

DEVELOPMENT OF ATTENUATION FORMULAS 


The attenuation of the shock wave in a shock tube is now 
derived. The flow is considered to consist of the ideal 
shock-tube flow plus small perturbations due to the boundary 
layer. The ideal flow is uniform in regions 1, 2, 3, and 4, 
and is denoted by the appropriate subscripts (i. e., p u p 2 , Ps, 
etc.). Perturbations at any point are denoted by A and an 
additional subscript indicating the point. Thus, the per- 
turbations at point b of figure 2 evaluated on the right side 
of the entropy discontinuity are designated by A p 2>b} A w 2i& , 
and so forth. The expansion wave of the ideal flow is 
assumed to have negligible thickness and to propagate into 
region 4 with the velocity u = — a 4 as indicated in figure 2. 


(The same assumption was used in the attenuation study of 
ref. 6.) 

Let point d of figure 2 represent an arbitrary point on the 
shock wave. The problem is to find the net pressure 
perturbation behind the shock (i. e., Ap 2>d ). This requires 
an integration of equation (4) along all characteristic lines 
which contribute to A p 2td . Because of the entropy dis- 
continuities, there are an infinite number of line segments 
along which the integration must be conducted (fig. 2). 
However, the major contributions to A p 2td can be shown to 
come from segments ab, cb, and bd. 
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If Ap£ d represents the incident pressure wave at point d, 
then equation (C7) gives 


Ap 2id =FAp£ d 


(FI) 


But, from equation (4), 


A ^= A *»+d?Tk) (F2) 

where the integration is conducted along the line bd. From 
equation (C3), 

&Pt.i,= CApi b +DApZ b (F3) 


Again, from equation (4), 


hplb ~ tpie+ ada-m X„ ® 2( * ,T)d * 


(F4) 


It can be shown that Ap£ a =0 in the present case. More- 
over, p 2 =p 3 . The expression for A p 2td can then be expressed 


as 

1 a 2 d Aflu f 1 C x 
F 2y 2 ^2 |_l+Af 2 Js 6 


%dH 


O 

1~M 2 



Y 32 F 

< 232 ( 1 +^ 3 ) 




In appendixes D and E it is shown that for wholly laminar or 
wholly turbulent boundary layers 



where L and n are independent of r and £. The values of L 
and n depend on whether the boundary layer behind a given 
wave is laminar or turbulent. Substituting equations (F5) 
into equation (5) and using equations (B3) yield 


(infer© 


n 2 A p 2 d 


P2 


2 t 2 FL 2 (u s la 2 y- 2n z 


1—^2 fi 

, 1— K/a 2 )T 2 

J 

m 2 J 


i+fC-Ts, 1 


-n 2 L 32 D (Af 2 ) 1+3( ”2~”3 ) 


02 ) 

M 2 / v 3 n 


1— n 3 a 32 

/at\*-** 

(l+M 3 +a43) r 2 “ 2 


\u s J 



y,n 2 — n 3 


(F6) 


The corresponding change in shock Mach number can be 
found from equation (C5) . Appropriate values for L and n 
are given in appendixes D and E. 

If n 2 =n 3 and CEAp 2 jAp 2fd is assumed not to vary with x, 
then equation (F6) indicates that Ap 2 jAp 2td ={% c lx) l ~ n * and 

CE^^=CE(^y ni (F7) 

A P2.a \xj 


Equation (F7) indicates CEAp 2iC /Ap 2 , d to be independent of 
x and is therefore consistent with the original assumption to 
this effect. Substitution of equation (F7) into equation 


(F6) gives accurate results even when n 2 9^n 3 , since the major 
contribution to A p 2yd comes from the first term on the right 
side of equation (F6). The term CEAp 2tC /Ap 2td represents 
the contribution to A p 2td of all the characteristic lines of 
figure 2 other than lines ab, be, and bd, and is referred to as 
the contribution of the reflected wave at point c in figure 6. 
It can be seen from figure 6 that this term is small and can 
probably be neglected in most cases. 

If it is necessary to consider the boundary layer behind 
the shock (or expansion) wave as partly laminar and partly 
turbulent, then equation (5) must be integrated accordingly. 


APPENDIX G 


IDEAL SHOCK-TUBE RELATIONS 


Shock-tube relations, assuming ideal flow, are presented 
herein for convenience. The formulas were obtained from 
reference 1. The notation Pi 2 ~PilP 2 , a 32 =a 3 ja 2 , and so 
forth, is again used. 


Define 


Then 


«= (yT - i)/Oy 1 ) 
0— (7 1)/2t 


B -s[ 1 -<»-«V« I rH IS]’' 


PU=(PuP2l) i/yi 


P 21 


l + «lff21 
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(Glc) 


T34 — a| 4 — Q?14#2l) 2 ^ 4 

rr 2 r^ 2 i(ai+^ 2 i)l 

T2i ~ an - L i+ssrj 


U 2 _ Pn—1 
0/1 Yl [@l(&lp2l~\~ 1)]^ 
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a* 7404 L J 

[ft(l+«tf«,i) J 

_ fti- 1 
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(Gle) 
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(Gli) 
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APPENDIX H 


REYNOLDS NUMBERS AND TRANSITION 


A Reynolds number characterizing the boundary-layer 
development in regions 2 and 3 of the shock tube is now 
developed. 


Region 2 (— <-L < l) 
° \ u s - u s t - ) 


In the notation of reference 2, a Reynolds number char- 
acterizing the boundary layer behind a shock is defined 
a _x_ ( u e~ u w) 2 
u. 


therein as 


Be=~ 


Using the transformations 


indicated by equations (Dl) and arbitrarily basing v on free- 
stream conditions, the Reynolds number for region 2 may 
be written as 


u 2 (u s r— Q u 2 /u s 

* 2 “ r 2 (1 -uju,) 


(HI) 


where £,r are the coordinates of a point in region 2 (following 
sketch) . 



Consider £,r to be a point on the characteristic line through 
x,t (above sketch). The Reynolds number for points along 
this line, as a function of £, is (using eq. (B2)) 


30 _u 2 (x— £) / u 2 ju s \ /T+M 2 — u s /a 2 \ 
2 ~ r 2 1 +M 2 ; 


(H2) 


If the transition Reynolds number is known, equation (H2) 
can be used to determine the values of £ at which transition 
occurs. The maximum Reynolds number in region 2 occurs 
at the contact surface (|=|*) and equals 

^•*=vW3 2 ( i+m -3^] (H3) 


Equation (H3) can be used as a general index as to whether 
the boundary layer (along the characteristic line) is primarily 
laminar or primarily turbulent. A plot of Re 2>b against M s 
for air in region 1 (with ri = 520° R) is given in figure 6. 

The transition Reynolds number for the boundary layer 
behind a shock wave has not yet been established. In the 
absence of more accurate information, the transition Reyn- 


olds number for incompressible flow over a semi-infinite 
flat plate [I?e=0(0.5X10 6 )] might be used to estimate the 
transition point behind a shock wave. That is, take 
#e 2 =0(0.5X10 6 ) as a rough estimate for the transition 
Reynolds number. This is probably a conservative estimate 
for the strong-shock cases, since the large amount of heat 
transfer to the shock-tube walls may have a very stabilizing 
effect on the boundary layer. 

Region 3^—1 <— ^ 

\ a A r aj 

For a point in region 3 (following sketch), the Reynolds 



number as defined in reference 2 becomes (using eqs. (D2) 
and v based on the free stream) 


p _ u 3 (a4T-\-& Us/ai 

v 3 1-f -u 3 (a\ 


(H4) 


For a point on the characteristic line influencing point x,t, 
the Reynolds number becomes 

D u 3 x u z ja A 1+A/8+<z 43 ["£ /l+A/ 2 — u s la 2 \a,r] 
v t l+uja* 1 +M 3 LA l+iH 3 + a Ju s J 


At the contact surface (£=£»), equation (H5) has the value 


(H6) 


r l'he criterion Re 3 — 0(0.5 X10 6 ) miglit be taken as a rough 
estimate for the transition Reynolds number behind the 
expansion wave. 

The boundary-layer characteristics presented in appendixes 
D and E indicate a discontinuity in the boundary -layer profile 
across the contact surface. (Note that 
This discontinuity does not actually occur, and the theoretical 
discontinuity thus represents a deficiency of the present 
method. 
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APPENDIX I 


GENERATION OF PRESSURE WAVES BY WALL SHEAR AND HEAT ADDITION 


The generation of pressure waves by body forces and by 
heat sources in one-dimensional flow is investigated, and 
the results are used to compare the method of the present 
report with that of reference 6. 

WAVE GENERATION BY BODY FORCES AND HEAT SOURCES 

Assume a uniform flow of pressure p, velocity u, and so 
forth, in a tube of constant cross section. This uniform 
flow is assumed to be slightly perturbed by weak body 
forces and heat sources. The equations of motion are 
(neglecting the possibility of mass sources, since this case is 
treated in eq. (1)) 


fbAu 

p Ur+ 

bAu\ 
' U bx / 

=/- 

bAp ^ 
bx 

dAp i D 

bAu . 

bAp 

— n v 

c>< +p 


bx 


bAs 


_ 2 


bt 

+ “ daT 

pT 



where f(x,t ) is the body force per unit volume acting in the 
+ ^-direction and q(x,t) is the heat addition per unit volume 
per unit time. The symbol As represents the entropy 
perturbation of a particle and is related to the pressure and 
density perturbations by 


(1) and (II) are linear, the solution for the case where mass 
sources are also present can be obtained by adding equation 

(2) to equation (13). Note that an elemental heat or mass 
source generates a symmetric pressure-wave pattern (posi- 
tive pressure waves propagating in the downstream and 
upstream directions) , while an elemental body force generates 
an antisymmetric pressure-wave pattern (positive pressure 
waves propagating in the downstream direction and negative 
pressure waves propagating in the upstream direction). 

COMPARISON WITH REFERENCE 6 

In effect, reference 6 uses equations (13) to find the shock 
attenuation in a shock tube. The value of g is obtained by 
averaging, across the tube cross section, the heat transfer at 
the wall and the viscous dissipation in the boundary layer. 
The value of / is found by averaging the wall shear across the 
tube cross section. Thus, if g w is the heat transferred into the 
boundary layer per unit wall area, H w is the net dissipation in 
the boundary layer per unit wall area, and r» is the shear per 
unit area exerted by the wall on the fluid (taken to be positive 
in the + ^-direction), then 

2=3 0„+ff«)=2 (a.+flj 

-i- 

t ~A Tw ~d T " 



As_ Ap A p 

c v ~ P P 


( 12 ) 


The perturbation at any point x,t can be shown to equal 

Ap_Ap + A p~ 

V 

Au__l / Ap + A p~\ 
a ~y \ p p ) 


As 

c 


C x 


■■ 


c n T 




Ap__l / Ap As\ 
P V \ p c v ) 


(13) 


where 

Ap + 


p 2pa(l-f-M)J_ 

A P~^ y 
p 2pa(l—M) 


L c,T ^ a 

rr«(«-HK) 'fa-K)’ 

Jx L e,T a 




The upper limit on the integral for A p~ is -)-<*> or — oo 
depending on whether M<^ 1 or AT> 1. Since equations 


Substituting equation (14) into equation (13) and integrating 
along appropriate characteristic lines should yield the same 
results for shock attenuation as those which were obtained in 
reference 6 from a somewhat different viewpoint. However, 
these results are not in agreement with the results obtained 
from equation (4) (the latter being the basis of the present 
report) . 

Equation (4) is based on v while equation (14) contains 
terms relating to the heat transfer at the wall, viscous dis- 
sipation, and wall shear. But, from boundary-layer theory, 
it can be shown that v is dependent on terms related to the 
heat transfer at the wall, viscous dissipation, and the velocity 
profile (see, e. g., eqs. (D3a) and (D3c)). The heat-transfer 
and dissipation terms in equations (14) and (4) play the same 
qualitative roles and, therefore, are not discussed further. 
However, the wall-shear term in equation (14) in some cases 
has a sign opposite to that of the velocity-profile term in 
equation (4) . Thus, these terms are not always in qualitative 
agreement, and this leads to discrepancies between the 
results obtained from equations (14) and (4). 

Consider, for example, the boundary layer in a shock tube 
for the case of negligible heat transfer and dissipation. For 
this case, v depends only on the velocity-profile term. Simi- 
larly, only the term is retained in equation (4). The signs 
of v and r w and of the resulting pressure perturbations in 
regions 2 and 3 of the shock tube are summarized in the 
following table: 
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Region 

Sign 

Perturbations based 
on v (eq (4)) 

Perturbations based 
on t U) (eq (14)) 

V 

A p + 

A p~ 

T w 

Ap + 

A p~ 

2 

3 

+ 

+ 

+ 

- 

- 

+ 

+ 


Thus, when heat transfer and dissipation are neglected, a 
perturbation solution based on v differs from a perturbation 
solution based on r w in regard to the sign of the A p~ waves in 
region 2 and the Ap + waves in region 3. 

It may be concluded that attenuation solutions based on 
equation (14) differ from those based on equation (4) in the 


following respects: (1) The integration of equation (14) along 
characteristic line ab tends to overestimate its contribution 
to shock attenuation, and (2) the integration of equations 
(14) along characteristic line be tends to underestimate its 
contribution to shock attenuation. For the limiting case 
M a -* 1, the integrations along ab and be are negligible, and 
equations (4) and (14) give the same results for shock 
attenuation. 

When the boundary layer is thin, it is obvious that the 
attenuation calculation should be based on equation (4). If 
the viscous shear affects the entire cross section (as for a long 
shock tube), there is no longer a core of potential flow and the 
solution should probably be based on equations (14) (or on 
the characteristics method of ref. 6 when nonlinearities be- 
come important) . The quantities 7«, q_w, and H w should then 
be based on unsteady pipe flow rather than on thin-boundary- 
layer theory. In practice, the boundary layer is generally 
sufficiently thin to permit the use of equation (4) . 
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17.025 

20.771 

24. 418 

28. 006 

62. 879 

97. 285 

131. 58 

165. 84 

200.08 

285. 63 

371. 15 

456. 65 

542. 16 

627. 65 

884. 13 

1739. 0 

.9 



12. 127 

15. 434 

18. 572 

21.630 

24.643 

53.951 

82. 880 

111. 72 

140. 53 

169. 32 

241.26 

313. 18 

385. 09 

456. 99 

528. 89 

744.57 

1463.5 

.95 



11.043 

13.600 

16. 050 

18. 447 

20. 813 

43. 926 

66. 777 

89.568 

112. 34 

135.09 

191. 95 

248.80 

305. 64 

362. 47 

419. 31 

589. 79 

1158. 1 

.975 



10.384 

12.488 

14. 525 

16. 529 

18. 512 

37. 979 

57. 267 

76. 514 

95. 744 

114. 97 

163.00 

211. 03 

259. 05 

307. 06 

355. 08 

499.12 

972. 39 

1.00 



0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 


CD 


ATTENUATION IN A SHOCK TUBE DUE TO UNSTEADY-BOUND ARY-LAYER ACTION 




